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Abstract. Spanning the whole functional space of cosmologies with any admissible 
DE state equations w(a) seems a need, in view of forthcoming observations, namely 
those aiming to provide a tomography of cosmic shear. In this paper I show that 
this duty can be eased and that a suitable use of results for constant-?/; cosmologies 
can be sufficient. More in detail, I "assign" here six cosmologies, aiming to span the 
space of state equations 10(a) = w a + w a (l — a), for w and w a values consistent 
with WMAP5 and WMAP7 releases and run N-body simulations to work out their 
non-linear fluctuation spectra at various redshifts z. Such spectra are then compared 
with those of suitable auxiliary models, characterized by constant w. For each z & 
different auxiliary model is needed. Spectral discrepancies between the assigned and 
the auxiliary models, up to k ~ 2-3/iMpc~ 1 , are shown to keep within 1 %. Quite in 
general, discrepancies are smaller at greater z and exhibit a specific trend across the 
w and w a plane. Besides of aiming at simplifying the evaluation of spectra for a wide 
range of models, this paper also outlines a specific danger for future studies of the DE 
state equation, as models fairly distant on the WQ—w a plane can be easily confused. 
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1. Introduction 

One of the main puzzles of cosmology is why a model as ACDM, with so many conceptual 
problems, fits data so nicely. It is then important that the fine tuning paradox of ACDM 
is eased, with no likelihood downgrade [TJ [2], if Dark Energy (DE) is a self-interacting 
scalar field (dDE cosmologies). 

Although several researchers privilege potentials allowing tracking solutions [3, 4], 
data on V(<p) can be recovered just by testing the evolution of the DE scale parameter, 
w(a). Here a = 1/(1 + 2) is the scale factor in the spatially flat metric ds 2 = 
c 2 dt 2 — a 2 (t)d£ 2 , with dl being the comoving spatial distance element. 

Using available data, the WMAP team [5j [6] tried to constrain the coefficients Wq 
and w a in the expression 

w(a) = w + (1 - a) w a (1) 

for the DE state parameter. However, the setting of the likelihood ellipse suggested, on 
the w -w a plane, has significantly changed from WMAP5 to WM AP7 release. The two 
ellipses are overlapped in Figure [TJ where we also indicate the boundary line w Q = —w a , 
beyond which the DE state equation should be rejected, unless further modified by other 
parameters at high z (otherwise, DE could become too dense, possibly modifying BBN 
and even Meszaros' effect). 

A further selection among models in such ellipses shall be provided, in a near future, 
by different observations (e.g. [HE]) and, in particular, by tomographic shear surveys 
(e.g. [9]), able to reconstruct the matter fluctuations spectra P(k, z) at various z's, with 
a precision approaching 1 % [10] . It is therefore important to provide a tool to ease the 
determination of model spectra with complex DE state equations; namely avoiding the 
need to explore the whole functional space of w(a). 

Quite in general, such spectra are to be obtained through N-body simulations; for 
cosmologies with a variable state parameter w(a), they have been performed since 2003 
(e.g. [II]. [12]. [13], [33]) and compared with w = const simulation outputs. Observables 
considered in these papers, however, only marginally included spectra. 

An important step forward was then due to Francis et al. [15] . They showed 
that suitable tuned constant-w models, at z = 0, closely approximate the spectra of 
cosmologies with a state parameter given by eq. ([I]). More precisely, the spectrum 
P(k, 0) of an assigned (A) model with state parameter w(a), can be approached by 
an auxiliary constant--?/; model (W) such that: (i) in A and W, Qb,m,tot, h and <r 8 are 
equal, (ii) the constant DE state parameter of W is tuned to yield equal comoving 
distances from the Last Scattering Band (LSB) and z = 0, for W and A. Then 
spectral discrepancies keep < 1 %, up to k ~ 2-3 h Mpc -1 (here symbols have their 
usual meaning). 

Discrepancies between A and W, at higher z, were also tested in [?], but they 
increase up to several percents, so that their precision may not be enough to exploit 
data. 
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Figure 1. Likelihood ellipses on the w -w a plane from WMAP5 and WMAP7 data 
releases. The black triangles indicate the A models considered in this work. State 
equations beyond the w a = —w a line should be modified at high- z. The dotted lines 
cross on the model best fitting recent data. 



The required precision was then obtained through a technique introduced by 
Casarini et al. [16] (Paper I hereafter) and tested for two specific cosmologies (see 
also [T7], for a further extension based on hydrodynamical simulation). Here I plan to 
test more cosmologies, so to sample the parameter space compatible with WMAP5 and 
WMAP7 data, also exploring the precision trend in its different regions. The setting of 
the six A models considered is shown in Figure [TJ The other parameters are consistent 
with both WMAP5 and WMAP7 data: matter density Q m = 0.274, Hubble parameter 
[100 km/s/Mpc] h = 0.7, fluctuation amplitude at 8/i _1 Mpc ct 8 = 0.81 and scalar 
spectral index n s = 0.96. 

The plan of the paper is as follows: in §2 is described the approach used in Paper I, 
§3 is devoted to describing our simulations and the techniques used to analyse them, in 
§4 are presented our results, and in §5 are discussed them. In Appendix A is reported 
the algebraic technique used to reproduce ellipses in Figure HJ 
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2. The spectral equivalence criterion 

Let me then first recall the technique presented in Paper I. At variance from [?], given 
an assigned model A, we introduced a specific auxiliary model W(z), for each z\ A and 
W(z), first of all, are required to share the values of w& )C , m = fib,c, m h 2 and erg, at such z . 

The former request is easily fulfilled; in fact, at any redshift the critical density p cr 
is defined through the value of the Hubble parameter H, being H 2 = (8nG/3)p cr . If we 
multiply both sides of this relation by Q m (or Qb, ^c) we have 

n m H 2 = (87iG/3)p m . (2) 

The r.h.s. of this equation, and then u m oc Q m H 2 (or Ub, u) c ), scale as a~ 3 , independently 
of the model. Accordingly, once A and W share Ub, c ,m at z = 0, it is so at any z: all 
W(z) models have equal Ub tC ,m- 

On the contrary, the evolution of erg depends on DE state equation. Its value at 
z = 0, as well as the value needed to normalize initial conditions, can be worked out 
only once we know the constant DE state parameters w(z) of the W(z) models. 

We come then to the most specific requirement, causing the dependence on z of 
the constant w's: that w is tuned so that W(z) and A have equal comoving distances 
between z and the LSB. 

The choice of H (the Hubble parameter at z — 0) is still unconstrained. Taking it 
however equal to H Q in A yields boxes with equal side L in both Mpc and /i _1 Mpc units. 
Notice that a simple-minded generalization of the criterion in [?], to high z, requires 
equal Q b ^ m (z) and, thence, H(z); this would create serious problems of sample variance 
and model comparison. 

3. Simulations 

Simulations performed for this work are meant to test the spectra of the M. models 
against the corresponding auxiliary W models up to z = 2. We compare simulations 
starting from realizations fixed by using an identical random seed. Initial conditions 
have been created, at z — 24, with the same procedure as in Paper I. They were then 
run by using the pkdgrav code [IE], modified to deal with any variable w(a) for Paper 
I. All models are run in a box with side L box = 256/i~ 1 Mpc, using N = 256 3 particles 
and a gravitational softening e = 25/i^ 1 kpc. 

Besides of the six models A, we have 4 auxiliary models W(z) for z = 0, 0.5, 1,2, 
which were run just down to the redshift where they are tested. Altogether, therefore, 
we run 30 model simulations. 

Model spectra are then worked out through a Fast Fourier Transform (FFT) of the 
matter density field. This last quantity is computed on a regular grid Nq x Nq x Nq 
(with Nq = 2048) from the particle distribution via a Cloud in Cell algorithm. 

Mass functions were also worked out for all models and found to be consistent with 
Sheth & Tormen predictions. 
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Figure 2. Each box in this Figure refers to a single A model (black triangle) and 
the related W(z) models (color polygones). Notice that, in most cases, distances 
between colored polygones are smaller than their distance from the black triangle. 
This outlines the possibility of a serious bias in data analysis, if tested by assuming 
constant w cosmologies. 



4. Results 

In what follows, A models will be ordered according to increasing values of the parameter 
w Q , from Al to A6. In Figure [2] I report the setting of *4's and related W-models on 
the w -w a plane, for z = 0, 0.5, 1, 2 . 

The Al cosmology was consistent with WMAP5 and is apparently outside the 2- 
er curve for WMAP7. It should be however reminded that the significant shift of the 
ellipses is unlikely due to the fresh CMB inputs, while omitting to impose the distance 
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Figure 3. Spectral discrepancies at various redshift. Each box refers to a single A 
model (same ordering as in previous Figure). Discrepancies are however greatest at 
z = and, at such redshift only, attain 2% in models with w -C — 1. The k value 
where such large discrepancy is reached, however, is ~ 4-5 h Mpc -1 , well above the 
fc-range where gas dynamics can be ignored, if aiming at ~ 1 % precision. Keeping 
within k ~ 4-5 h Mpc -1 , the top discrepancy is always < 1 % . 

prior [51 16] surely had an impact on it. 

Figure [2] indicates that the distance between W models, lying on the w a = line, 
by definition, is mostly smaller than their distance from A. Distances however scale 
with w a and are smaller for the central w Q values. I shall return on this point in the 
next Section. 

Figure [3] then shows spectral discrepancies. Within k = 3 h Mpc -1 the maximum 
discrepancy is reached for Al and A2 at z = 0, attaining 1.4 and 1.6%, respectively. 
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Owing to [T7], however, this is a scale where baryon dynamics affects spectra already 
more than 1-2%. In Paper I we actually took k = 3/iMpc" 1 as a limit; in this 
scale range, however, the discrepancy is rapidly bursting and, should we keep within 
k = 2 /iMpc -1 , no spectral discrepancy exceeds 1 %. 

It is however clear that the other models exhibit a nicer behavior. Even the A6 
model, whose distance from ACDM is similar to Al, exhibits discrepancies within 0.8 %, 
up to k = 10 . Other models are even nicer, exhibiting discrepancies in the range of the 
permil. 

Discrepancies decrease with increasing z and, at z = 0.5 are mostly in the permil 
range, for all models. 

5. Discussion 

This finding perfectly fits the reason why spectral similarities are expected. When using 
the conformal time r, the background cosmic metric reads ds 2 = a 2 (r)(dT 2 — d£ 2 ). Equal 
comoving distance means then that equal conformal times have elapsed. Most events, on 
a cosmological scale, are indeed scheduled in accordance with the conformal time r ; the 
ordinary time t, instead, sets a correct "timing" in virialized environments, where local 
minkowskian reference frames no longer feel the scale factor evolution. Accordingly, 
it makes sense to set aside models with equal "conformal age", expecting increasing 
discrepancies when time elapses and, however, on scales virialized since longer. Up to 
k ~ 2-3 h Mpc" 1 we are inspecting scales >~ L = 27r/k ~ 2 /i _1 Mpc; still at z — 0, 
there are quite a few virialized systems on such scales, which however correspond to 
density peaks more and more above average, as we go towards higher z . 

A further comment is deserved by the proximity of auxiliary models on the w -w a 
plane. In quite a few CclSGS, clS for the models A3 and A4, they almost overlap. This 
envisages a danger, in future observational analysis; it is reasonable to expect that a 
first data test is carried by assuming constant w. Let us suppose, for instance, that the 
real cosmology is close to the model AA. It is not unlikely that all the models indicated 
by colored polygones are then compatible with a single w value, e.g. -1.18±0.10 . 

As a matter of fact, starting from the setting of each triangle, we could draw a bunch 
of curves, indicating the loci of equal r — r rec (difference between conformal present and 
recombination times), when w a varies. When these curves diverge fast enough, there 
is a realistic possibility that they cross the w a = line on reasonably distant sites. 
Otherwise, the risk of spurious constant-//; detection is a serious danger. One should 
also take into account that the very assumption of a polynomial w(a) is a simplifying 
ansatz, and that there can be cosmologies behaving even more dangerously, e.g. some 
cosmology arising from a tracking potential. 

In this context, the very discrepancies detected above k ~ 2-3 h Mpc -1 are welcome. 
These scales still need to be tested through hydro simulations, but a reliable pattern for 
data analysis could actually start from the assumption of constant--?/;, so individuating 
a bunch of curves, characterized by constant r Q — r rec , which will be the models among 
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which one will discriminate through higher-A; spectral discrepancies. 
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Appendix A. Reproduction of the likelihood ellipses with the Bezier curves 

The algebraic technique used to reproduce the ellipses in Figure Q] is named after Bezier 
and is largely used in vector graphics to model smooth curves which can be scaled 
indefinitely, without any bound, by the limits of rasterized images. In the PostScript 
files in the WMAP A-site, I found the coefficients for the cubic Bezier curves: 

B(u) = (1 - u) 3 P + 3(1 - ufuP l + 3(1 - u)u 2 P 2 + u 3 P 3 u E [0, l](A.l) 

yielding 1- and 2-cr contours. Here the vector B, running on the wo-w a plane, describes 
a curve fixed by the positions of the points (k=0,...,3), when u varies from to 
1. Further details on this technique can be found in a previous paper [19], where the 
coordinates of Bezier points for WMAP5 ellipses are also reported. Here below I report 
the coordinates of the Bezier points to draw WMAP7 ellipses. 
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Table Al. Points defining the 4 cubic Bezier expressions yielding the 1-a and 2-cr 
curves (upper and lower table, respectively). 
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